Sterile neutrino decay and the LSND experiment 
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Abstract. We propose a new explanation of the intriguing LSND evidence for electron 
antineutrino appearance in terms of heavy (mostly sterile) neutrino decay via a coupling with 
a light scalar and light (mostly active) neutrinos. We perform a fit to the LSND data, as well 
as all relevant null-result experiments, taking into account the distortion of the spectrum due 
0^ ■ to decay. By requiring a coupling g ~ 10~^, a heavy neutrino mass ^ 100 keV and a mixing 

with muon neutrinos |(7^i4p ~ 10~^, we show that this model explains all existing data evading 
constraints that disfavor standard (3-1-1) neutrino models. 



The confirmation of tlie solar ^ and atmospheric neutrino oscillations by reactor [Hj and 
. accelerator ^ experiments and the bounds from laboratory neutrino mass measurements jH] have 

I set the "standard" neutrino picture in terms of three light active neutrinos with two mass square 

' differences. However, the observed excess of i>e in the LSND experiment [S] was interpreted as 

evidence of z?^ — > transitions with a mass square splitting, Am'j^gj^j^ ^ 1 _ 10 eV^, that 
cannot be accommodated by the other two, Am^^^ ~ 8 x 10~^ eV^ and Am^^^ ~ 2.5 x 10^^ eV^. 
The most natural solution would be to introduce a fourth neutrino |7| with the appropriate 
mass to give the hinted mass square difference. However, this does not turn out to lead to a 
satisfactory description of all data in terms of neutrino oscillations |H1 because of tight constraints 
. from solar [J, atmospheric PJ, and null-result short-baseline (SBL) experiments [SJ. In view of 

this, during the last years, other alternative explanations have been proposed |lUllllj . Unlike in 
Ref ^01 where decay was invoked to evade SBL results, here we propose (see Ref. ^2] for details) 
a new explanation of LSND in terms of heavy (mostly sterile) neutrino decay. This new massive 
' state, n4, must have a small mixing with the muon neutrino, and the signal is explained 

by its decay into a combination of light neutrinos, being predominantly of electron type. 

A natural way to introduce neutrino decay is by means of a term in the Lagrangian, which 
couples neutrinos to a light scalar. We assume for the scalar mass mi^2,3 ^ n^(i> *S ^4, where 
"^1,2,3) nT-fj) and 7714 are the masses of three light neutrinos, of the scalar and of the new massive 
state, respectively. In this way the three light neutrinos are stable. Hence, the terms in the 
Lagrangian which provide the decay are given, in the mass basis, by 
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where / = 1, 2, 3 and, in general, the coupling matrix g/^i is complex. 

In the case of Majorana particles, neutrinos are identical to antineutrinos. Weak interactions 
couple to left-handed (chiral) neutrinos and right-handed (chiral) antineutrinos while the 
propagation states are those of definite helicity. In the relativistic case, one can identify neutrinos 
and antineutrinos with helicity states up to terms of order m/Ey. It follows from Eq. that 
not only the decay n ^ v + cf) can occur, but also n ^ v + 4>\s possible ^j. Hence, the expected 
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Figure 1. L/E^ spectrum for LSND for decay 
and oscillations compared with the data given in 
Fig. 24 of Ref. 0. The hatched histo gram shows 
the contribution of the line from the primary pion 
decay in the decay scenario. From Ref. |12|. 




Figure 2. Allowed regions for LSND + KARMEN 
(solid) and SBL disappearance+atmospheric neutrino 
experiments (dashed) at 99% CL, and their combi- 
nation (shaded regions) at 90% and 99% CL. From 

Ref. na. 



number of V(, events with neutrino energy in the interval [Ep^, Ep^ + dEp^] is given by |12j 
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where (l)^{Ep^) is the muon spectrum, 4>o = J dEp^(j)^{Ep^) and dPyr^^ys^{Ey^) / dE^^^ is the 
differential probability that a neutrino of flavor a with energy Ejj^ is converted into an neutrino 



of flavor (3 with energy in the interval , Ey^ 



+ dE, 



The indexes r, s take the value 



('+') for neutrinos (antineutrinos). The detection cross section is given by a{Ep^) and C is an 
overall constant containing the number of target particles, efficiencies and geometrical factors. 

Unlike other explanations of the LSND result based on sterile neutrinos, our model does not 
require mixing of the electron neutrino with the heavy mass state, so for simplicity we assume 
UeA = 0. Hence the differential probability for v'' i/| conversion can be written as |12j 
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where r4 is the decay rate of as follows from Eq. jl3j . 

Our analysis of LSND data (Fig.^, which does not use the much less significant decay in flight 
(DIF) sample, gives a best fit value of Xmin ~ 5.6/9 dof for oscillations and Xmin ~ 10-8/9 dof 
for decay. The reason for the slightly worse fit for decay is the spectral distortion implied by the 
energy distribution of the decay products. Nevertheless, the overall goodness of fit (GOF) is still 
acceptable, 29%, giving as best fit values |?7/j4p = 0.016 and gm^ = 3.4eV, where = J2i l5'4iP- 

As is well known, (3+1) neutrino models are disfavored due to the tension between LSND 
and SBL experiments, which have not observed either appearance or dissapearance so far 
We have also studied their compatibility in the decay scenario. Under the assumption UeA, = 0, 
no Ue disappearance is expected in SBL reactor experiments. However, since the explanation of 
the LSND signal within this scenario requires mixing of with the heavy state, one expects 
some effect in disappearance experiments. The main constraints on the allowed LSND- 
KARMEN region come from the CDHS experiment [0] and from atmospheric neutrino data |14j . 
Nevertheless, the different expression for the z/^ survival probability in the decay scenario 
provides a rather poor sensitivity to which comes via |[/^4|^, unlike the case of oscillations, 
that it comes via |[/^4p. The results are shown in Fig. |2l where appearance and disappearance 
experiments are shown to be in perfect agreement, with large overlap of the allowed regions. 



Mixing between active and sterile neutrinos, and couplings between active neutrinos and 
a light scalar, have been extensively studied, both in laboratory experiments and, for their 
implications in the evolution of the early Universe and of supernovae. The decay of pions and 
kaons has been used to set bounds on the mixing of heavy neutrinos for masses above MeV |15|| . 
In our model, the most stringent bounds on the coupling from laboratory experiments, coming 
from light scalar emission in pion and kaon decays, are of order (f' < few x 10~^ 16]. On the 
other hand, a coupling g ~ 10~^ would imply that and (j) are strongly coupled to active 
neutrinos inside the supernova, and hence, they are trapped within the neutrinosphere, avoiding 
any energy loss due to particles escaping from the core. This value of the coupling would also 
imply that as soon as are produced in the Early Universe, they are thermalized along with 
the scalar, leading to 1.57 extra number of relativistic degrees of freedom at BBN. 

In summary, we have presented an explanation of the LSND evidence for z^e appearance based 
on the decay of a heavy (mostly sterile) neutrino into a light scalar particle and light neutrinos. 
It is needed |C//x4p ~ 10~^ and gru/^ ~ 1 eV, and hence for g ~ 10~^, ~ 100 keV. In addition, 
the decay model also predicts a detectable signal in the upcoming MiniBooNE experiment \Vf\. 
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